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Abstract. The moduli space of Frobenius manifolds carries a natural involutive symmetry, 
and a distinguished class - so-called modular Frobenius manifolds - lie at the fixed points 
of this symmetry. In this paper a classification of semi-simple modular Frobenius manifolds 
which are polynomial in all but one of the variables is begun, and completed for three 
and four dimensional manifolds. The resulting examples may also be obtained from higher 
dimensional manifolds by a process of folding. The relationship of these results with orbifold 
quantum cohomology is also discussed. 



1. Introduction 

The simplest family of Frobenius manifolds are those for which the prepotential - a solu- 
tion of the underlying WDVV-equations - is a polynomial in the flat variables and the Saito 
construction [20] provides a construction of these polynomial solutions on the orbit space 
C'^/VF , starting from any finite Coxeter group W . It was conjectured by Dubrovin [7J and 
later proved by Hertling [12j that all polynomial solutions arise in this way, and hence polyno- 
mial Frobenius manifolds have now been completely classified. Beyond polynomial solutions 
there are trigonometric and elliptic solutions and variations of the original Saito construction, 
corresponding to extended affine Weyl and Jacobi groups, which provide large classes of ex- 
amples. One may conjecture that these constructions provide all such solutions within their 
classes, but currently, unlike in the polynomial case, such conjectures remain unproved. 

In [15] the authors introduced the idea of a modular Frobenius manifold. These sit at the 
fixed points of a natural involutive symmetry that exists on the moduli space of Frobenius 
manifolds and hence inherit very special properties, basically an invariance under the modular 
group. 

The aim of this paper is to study such modular Frobenius manifolds further and in partic- 
ular to begin a low-dimensional classification of polynomial examples, where this means that 
the prepotential is polynomial in the variables , . . . and transcendental in the last 

variable (throughout this paper N will denote the dimension of the Frobenius manifold). 
A number of explicit examples have appeared in the literature before: 

• the Ai Jacobi group example [7]; 

• the dJ^'^^ solution [21]; 

• the simple elliptic singularities Eej^s [131 ESI ED [231 ED] • 

What is perhaps surprising is that the simplest family of examples - those based on the Ai\f-2- 
Jacobi group - do not fall into this class, apart from the simplest A'^ = 3 case. In general these 
are polynomial in the variables , . . . t^~'^ but have rational dependence in the variable t^~^ 
[3]. Thus polynomial modular Frobenius manifolds seem to be extremely special. 



Date: l?**" Oct. 2011. 

1991 Mathematics Subject Classification. 53D45,14H70,11F03. 
Key words and phrases. Frobenius manifolds, modular functions. 

1 



2 



EWAN K. MORRISON, IAN A.B. STRACHAN 



Some other solutions may be found in [10] and in the work of Bertola [3] , and an aim of this 
paper is to provide a systematic analysis of the solution space of polynomial modular manifolds 
rather than just isolated examples. In fact a close examination of Satake's construction [23] 
(itself based on the Saito construction [21]) for the so-called codimension one cases shows 
that the prepotential will be polynomial and modular. Since there are only a finite number 
of codimension one examples, a plausible conjecture would be that semi-simple polynomial 
modular manifolds are actually finite in number. 

2. Modular Frobenius manifolds 
We recall the basic definitions and symmetries of a Frobenius manifold [7]. 

2.1. Frobenius Manifolds and the WDVV Equations. Frobenius manifolds were in- 
troduced as a way to give a geometric understanding to solutions of the Witten-Dijkraaf- 
Verlinde-Verlinde (WDVV) equations, 

d^F d^F d^F d^F 

^ ^ df^dtf^dt^ ' df'dt^dt^ dt^dtf^dt^ ' ^t^'^t'^^t'r 

for some quasihomogeneous function F[t). Throughout this paper ■q'^^ will be defined via 
fl'^^Vpn = where 

is constant and non-degenerate. We recall briefly how to establish the correspondence between 
Frobenius manifolds and solutions of WDVV. 

Definition 1. The triple (A, o,r/) is a Frobenius algebra if: 

1. (A, o) is a commutative associative algebra over C with unity e; 

2. The bilinear pairing rj and multiplication o satisfy the following Frobenius condition 

r]{X oY,Z) = r]{X, Y o Z) , X,Y,Z eA. 
With this one may define a Frobenius manifold. 

Definition 2. Let M be a smooth manifold. M. is called a Frobenius manifold if each tangent 
space TfM. is equipped with the structure of a Frobenius algebra varying smoothly with t £ A4, 
and further 

1 . The invariant inner product r] defines a flat metric on Ad . 

2. The unity vector field is covariantly constant with respect to the Levi-Civita connection 
for T], 

(3) Ve = 0. 

3. Let 

(4) ciX,Y,Z) ■.= r]{X oY,Z), X,Y,ZGTtM, 

then the (0, 4) tensor Vwc{X, Y, Z) is totally symmetric. 

4. There exists a vector field E G T{TAi) such that WE = and 

(5) CeV = (2 - d)r], Ceo = o, Cec = -e. 
E is called the Euler vector field. 



POLYNOMIAL MODULAR FROBENIUS MANIFOLDS 



3 



Condition 1 implies there exists a choice of coordinates {t^,...,t^) such that the Gram 
matrix ija^ = {da^dp) is constant. Furthermore, this may be done in such a way that e = 
di. In such a coordinate system, partial and covariant derivatives coincide, and condition 3 
becomes Cap-y^K = Capn^'y- Successive applications of the Poincare lemma then implies local 
existence of a function F{t) called the free energy of the Frobenius manifold such that 

- ^^-^(^) 

Since r]{X, Y) = ri{e o X,Y) = c{e, X,Y), we have 

(7) ValS = Clap- 

Defining (t/q,^)~^ = the components of o are given by c^^ = rj^^Ci^p^. Associativity of 
o is then equivalent to ([TJ. We assume from now on that VE is diagonalizable and that 
rjij = 5iJ^j^N+\ ■ Condition 4 leads to the requirement that F is quasihomogeneous, that is, 

(8) CeF = E(F) = dp • F = {3 — d)F modulo quadratic terms, 

and, on using the freedom in the definition of flat coordinates, the Euler vector field may be 
taken to be 

(9) E{t) = Y,{'^-^-^i^)e^^+ 

a {a:qc,=l} 

We will, for the most part, restrict ourselves in this paper to those Frobenius manifolds with 
the property that r'^ = , cr = 1 , . . . ,N. So 

(10) E = ^(1 - ^ - ^„)ra„ , 

a 

(11) CeF = (3-d)F. 

This excludes those examples coming from quantum cohomology and the extended affine Weyl 
orbit spaces |9]. 

To explain the origins of modular Frobenius manifolds it is first necessary to consider 
symmetries on the moduli space of Frobenius manifolds [7]. 

Definition 3. (a) Legendre-type transformations S'^; This is defined by 

ta = dtcdf^Fit) , (note : = ija/si'^) 
d^F _ d'^F 
WdtP ~ dt^dtP ' 

(b) Inversion Symmetry /: This is defined by 
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In this paper we will concentrate on inversion symmetries (also called type II symmetries). 
For these, the structure constants of the inverted Frobenius manifold are related to those of 
the original by 

(13) c„,, = it )- ^^^c.,.. 

Under the assumption that = for i = 1 , . . . ,N the corresponding Euler vector field of 
the inverted Frobenius manifold F has the form 

(14) E(i) = T.(^-^-f^a)i''d^, 

a 

where 

(15) d = 2- d, /ii = /iAT - 1, fiN = fJ.i + 1, fii= fJ-i, for 1,N. 

To motivate what follows, consider the following example (taken from [7J Appendix C). 
Example 4. Suppose N = 3 and let the Euler vector field and prepotential take the form 

(16) E = t'^ + it^A, F(^^^^^3) ^ i(,i)2,3 + i,i(,2)2 _ i^^^^sy 

The WDVV equations are equivalent to the nonlinear differential equation (known as the 
Chazy equation) 

(17) 7"' = 677" - 9(7')' 

for the function 7(t^) (here ' denotes differentiation with respect to the variable t^). The 
equation has an SL{2,C) invariance, mapping solutions to solutions, 

(18) t^ = ^(t^\ = (ct^ + dfj(t^\ ^2c(ct^ + d), ad-bc=l. 

ct-^ + d 

A simple calculation shows that under the inversion symmetry 

F{i) = l{t^)H^ + -i\ff-^^iJ-^(^)-i 

^ ^ 2^ ' ^ 2 ^ ' 16 \(t3)2 ^\t3 J t3 

where, on using the symmetry given by U8\} . 7 is also a solution of the Chazy equation. Thus 
under inversion symmetry one has a weak modular symmetry: the functional form of the 
prepotential is preserved, but with the function 7 being replaced by a new solution of the same 
equation connected by a special case of the symmetry [7gj) . 

One simple explicit solution of the Chazy equation is given by 

7(t3) = ^i?2(t^), 

where E2 is the second Eisenstein series. With this explicit solution one has a stronger 
symmetry. Since E2 has the modularity properties 

E2{t + 1)=E2{t), E2{-^=T^E2{t) + ^. 
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then under inversion symmetry one has F{t) = F{t) , i.e. both the functional form of the 
prepotential and the specific solution of the Chazy equation are preserved. 

The above example provides a motivation for studying those Frobenius manifolds with similar 
properties. 

2.2. Modular Frobenius Manifolds. 

Definition 5. Let M and M be two Frobenius manifolds with fixed coordinates t and t , with 
Euler vector fields E and E given by and , and connected via the inversion symmetry 
given by Definitions^ part (b). 

The manifold M is defined to be a modular Frobenius manifold if 

(19) F{i) = F{i) 
and 

E{i) = E{i). 

Thus instead of transforming the original prepotential into a different prepotential, the in- 
version symmetry maps the prepotential to itself, as in Example [H and hence one may think 
of modular Frobenius manifolds as lying at the fixed point of the involutive symmetry / . A 
comparison of the two Euler fields, using (jl4p and (jl5p gives the following constraints: d = d 
and hence d = 1 . Thus we have the following necessary conditions for a Frobenius manifold 
to be modular: 

(20) E = J^(--^,)r9„, 

a=l 

(21) CeF = 2F . 

Since the variables t^ and t^ behave differently to the remaining variables it is useful to 
use a different notation, namely: 

M = t\ = t\i = 2,... ,N -1, T = t^ , 

and z = (z^ , • • • , z^^^) . The two notations will be used interchangeably in what follows. 

A simple calculation then yields the modularity properties of the non-cubic part of the 
prepotential. 

Proposition 6. The prepotential 

(22) i7= l^V-l^x(z,z) + /(z,T) 

(where (x, y) = X^fj=2 ^ija^'y-' j satisfies (jl9|) under the action of I if and only if 

(23) /(^,-i)=l/(z,r)-^(z,z)^ 

Remark. There exists another choice of coordinates on a semi-simple Frobenius manifold, 
namely one that simplifies its algebraic structure. In these coordinates the multiplication is 
trivial: 

d ^ _ ^ ^ 
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It turns out that these coordinates are the roots of 

(24) det(g"''(t)-Ar/"^) =0, 

where g is the intersection form of the Frobenius manifold. The roots of the expression (j24|) 
are invariant under the symmetry /, and so the canonical coordinates are preserved up to 
a re-ordering [7]. This also provides a way to check whether the algebras are nilpotent or 
semi-simple: one just has to calculate the discriminant of this polynomial to check whether 
it has multiple roots (the nilpotent case) or not (the semi-simple case). 



3. Modularity and quasi-homogeneity 
We will consider polynomial prepotentials of the form (j22|) with 

(25) / = E I 11(^^4 5.(i^), 

aeL I i=2 ) 

SO z = {t^ ,T = t'^ and L := Z^J"^ . We do not assume any general properties of the 

functions Qa at this stage - the aim is to obtain the differential equations that these functions 
must satisfy in order for the WDVV equations to hold. Boundary or other conditions (for 
example, solutions being analytic at infinity) will then place constraints on these functions, but 
a priori no such constraints will be imposed. The following definition, while not immediately 
obvious, will play an important role in what follows. 

Definition 7. A pivot-point a € L is a lattice point for which 

coeff„ I J] ?7i,tV 1 y^O, 
where coeSa{p) is the coefficient o/ J3^2^(i*)°' of the polynomial p . 

Proposition 8. Assume that f defined above satisfies the modularity condition (I23p . Then: 
• If a is not a pivot-point then: 

5J--l=r(S«0-V(r); 



If a is a pivot-point then: 




9a {-^^ = r'^gair) - iTcoeffo, 

The proof is by direct computation and will be omitted. 

Having determined the modularity properties of the functions we now apply quasi-homogeneity. 
This will determine the possible terms that can appear in the ansatz (|25p . Recall that the 
Euler field takes the form 

E=y d,f^ 

1=1 



POLYNOMIAL MODULAR FROBENIUS MANIFOLDS 



7 



where di + d^-^-i-i = 1 . We now assume further that the di are positive rational numbers. 
Applying the quasihomogeneity condition E{F) = 2F (recall modular Frobenius manifold 
must have d = 1) implies the following constraint on the Oj : 

N-l 

(26) ditti = 2 . 

i=2 

Thus given the di we arrive at a special Diophantine equation whose solutions determines the 
possible monomials in (|25p (this special type of Diophantine equation is known in number 
theory Frobenius equation). 

Pivot-points - which are defined without reference to the Euler vector field - play an im- 
portant role in the construction of solutions of this equation. 

Lemma 9. Pivot-points automatically satisfy the Frobenius equation (|26p . 

Proof. On expanding ^X]fj^2 Vijt^'t^^ one can obtain the form of the pivot-points, namely 
(and here only non-zero elements are shown, all other elements are zero): 

(i) a = (. . . , 2, . . . , 2, . . .) with 2 in the i and N + 1 — i positions; 

(ii) a = (...,!,... ..,!,...) with 1 in the i , j and N + 1 — i , N + 1— j positions 

If is odd one obtains two further pivot-points: 

(iii) a = (. . . , 4, . . .) with 4 in the middle position; 

(iv) a = (. . . , 1, . . . , 2, . . . , 1, . . .) with 1 in the i and N + l — i positions and 2 in the middle 
position. 

Since di + dN+i-i = 1 (and hence (i(7v+i)/2 = 1/2 if is odd) the result follows. 
Given the above forms of the pivot-points one can easily count the number of such points: if 
is even the number of pivot-points is N{N — 2) /8 and if N is odd the number of pivot-points 
is (iV + l)(iV-l)/8. 

□ 

One can now give the geometric reasoning behind the name 'pivot-point'. The Frobenius 
equation ()26p defines a hyperplane 11 in a-space and we wish to find integer solutions, i.e. 
the points in 11 fi L . Since pivot-points are independent of the di as the di vary the plane 
n 'pivots' around these points. In the simplest non-trivial example, when = 4 , we have 
(1 — d3)a2 + d^a^ = 2 and a single pivot-point (2,2) . Thus as d^, varies the line rotates, or 
pivots, about this point: 

as 





(1 


-ds) 


a2 + c 


has = 


2 
















(2,2) 






d3 = 


2/3 






d3 = 


l/2\ 
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It is also obvious geometrically from the direction of the normal vector to 11 that the 
number of solutions to the Frobenius equation is finite and hence ()25p is a polynomial in the 
variables {^^ , . . . ,t^~^} . 

The number of independent pivot terms can be reduced further; in fact to one. Let 7(t) 
be any function with the transformation property 

J = t^i{t) - \t 

and let 

/(z, r) = 7(r)(z, zf + giz, r) . 
Then equation ([23|) implies that 

Thus one can obtain a refinement of Proposition [6) 

Proposition 10. The prepotential of a modular Frobenius manifold takes the form 

F = ^u^T - ^u{z, z) + 7(r)(z, zf + g{z, r) 

where 

T (-7) = ^'^W - ' 9 = r) . 

Moreover, if 

cfGLnn I 1=2 ) 

then 

ga (-i)=r(S°^)-V(r). 

Thus modularity and quasi-homogeneity determine the modularity properties of the func- 
tions 7 and ga together with the form of the monomial coefficients of these functions. To fix 
these functions one must now solve the WDVV equations (and by construction any solution 
will give a Frobenius manifold). 

Remark. Once the Euler vector field is fixed, the solutions of the Frobenius equation (j26|) 
determine the form of the prepotential. Imposing the WDVV equations then gives systems 
of over-determined non- linear ordinary differential equations (in the variable t^). It is these 
systems, determined by the specified di (subject, of course, to the constraint di + d]\f+i-i = 1) 
that are studied in this paper. These equations will possess similar properties to the Chazy 
equation in Example [H their solutions will have an SL(2, C) symmetry. To get a modular 
Frobenius manifold one then looks for solutions which have an invariance under SL{2,Z). In 
addition, if one then imposes the condition of semi-simplicity on the resulting multiplication 
one finds that this then leads to very strong constraints on the possible di . 

Another functional class in which to look for modular manifolds could be obtained by 
weakening the polynomial ansatz and replacing it with a rational ansatz (e.g rational in the 
variable t^~^ and polynomial in the variables t^ , . . . - this would then include the 

An -2 and -Bjv-2 examples of Bertola [3] . These examples, while rational, are constrained by 
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the condition that the functions gaiT) Sire never of negative weight, i.e. X^Oj > 2. Such a 
generahzation will not be pursued here. 

4. Solutions of the WDVV equations 

Proposition [10] gives the forms of the prepotential of a modular Frobenius manifold: the 
precise form of the modular functions that appear cannot be obtained from this analysis - 
this is fixed by imposing the WDVV equations of associativity. These equations decompose 
into three classes which depend on the number of r-derivatives. Consider the obstructions to 
associativity: 

A[X, Y,Z] = {X oY)oZ - X o{Y oZ). 
Since in the case being considered we have a unity element these simplify further: if any of the 
vector fields are equal to the unity field then A vanishes identically. Since the vector field dr 
is special (for example, it behaves differently to the other variables under modularity trans- 
formation), we decompose these equations further, taking the inner product with arbitrary 
vector fields to obtain scalar-valued equations. 

Proposition 11. The WDVV equations for a multiplication with unity field are equivalent 
to the vanishing of the following functions: 

A(^)(u, v) = r]{dr o dT-,u.o\^) — l](dr ou,dr ov) , 

A(^)(u, v,w) = 77(9,- o u, V o w) — 77(5,- o w, u o v) , 

A*^'^^(u, v,w,x) = 77(u o V, w o x) — ?7(u o X, V o w) 
for a// u , V , w , X € span{9ji i = 2 , . . . ,A^ — 1}. 
In terms of coordinate vector fields these conditions are: 

^ijk — {Vjk C-rri 'Hij C-rrk} ~\~ V^'^ {Cfip Cjkq Cjjg} ; 

(3) 

where r^ij = -{di,dj) . 

By judicious choice of variables one can obtain first order evolution equations for the 
unknown fields. However these systems are over determined, so one has to analyze various 
cases and subcases very carefully. 

4.1. = 3 solutions. In this simplest case (see Example H]) there is no freedom: the Euler 
vector field is fixed by the fact that d = 1 and hence 

^ id 1 2d 

The Frobenius equation (|26p only has one solution and hence modularity and quasi-homogeneity 
imply a polynomial prepotential of the forrcQ 

F = + Ihtl - ^47(^3) . 



^In examples only we lower indices for notational convenience, i.e. in expressions for the prepotentials in 
specific examples, U = . 
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a 


Non-pivot solutions (02,03): 


1 

2 




1 

3 


(3,0), (1,4), (0,6); 




(l,n+l),(0,2n); 


2n+l>(^>2) 


(0,2n + l); 


arbitrary 


none . 



Table 1. Non-pivot solutions of the Frobenius equation ()26p for = 4. 



The WDVV equations then imply that 7 must satisfy the third order equation 
(27) 7"' -677" +9(7')^ = 0. 

This is nothing more than the Chazy equation, whose modularity properties are well known 
(see [7]). This equation falls into Chazy 's class (see appendix) and its solutions may be written 
in terms of the Schwarzian triangle function S[^, g,0, t] . 

4.2. = 4 solutions. In this case there is a 1-parameter family of possible Euler fields, 
namely 

,id , ,^ ,.^d 



Without loss of generality we may assume decreasing degrees and hence a < ^ . A detailed 
analysis of the Frobenius equation ()26p gives the solutions summarized in Table 1. 

In all cases except o" G {| , |} the resulting solution of the WDVV equations gives a nilpo- 
tent algebra - one obtains Euler-type differential equations which may be solved explicitly, 
and the solution checked to see if the family of algebras it defines is semi-simple or not. The 
two remaining cases require a more detailed analysis. 

4.2.1. N = A ,a = ^ . In this case we take 

F = U\u + tit2h - \{t2hfi{U) + {4 gi{U) + ^2*1 52(^4) + tlt^ gz{U) + 4 54(^4)} • 

The analysis of the WDVV equations splits into two subcases (there are other subcases which 
appear in the analysis, but these yield non-semi-simple solutions): 

Case I: 52 = , 53 = 

One finds that 54 = ngi where /i is a constant and the evolution of gi and 7 are given by the 
equations: 



(28) 



7 
9'{ 



-7 - 288/i(7i, 
37 5i - 3517' • 



On eliminating gi one obtains a third-order scalar equation 
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(29) ■i = l<t^ + 877-10f. 

2 7 — 7 

Here the independent variable has been rescaled, t = , and hence 7' = ^7 etc.. This falls 
within Bureau's class (see appendix) and its solutions are given in terms of the Schwarzian 
triangle function y{t) = S[^, j,0,t] , namely 

7(t) = 




9iit) = . 

l92lJ,2 y2{y 

An alternative way to solve (|29p (following Satake [21] and section [6]) is to express the 
solutions in terms of solutions to the Halphen system 

UJi = -(^2^3 + '^l('^2 + W3) , 

L02 = -cJa^Ji + 072(^3 + wi) , 

Cjs = -U)iU2 + OJ-i{ljJi + OJ2) . 



One may show that the solution is: 

l{t) = ^(Wl +2CJ2 +^^3) , 
gi{t) = —^{u)i - 2uj2 + ws) . 

The required modularity properties of 7 and gi then follow automatically from the known 
modularity properties of the solution to the Halphen system. 

CaselL 52 /0,53 y^O 

One finds that 32 = , 93 = Kgi where K is a constant and the evolution of gi , g^ and 7 
are given by the equations: 

7' = ^7' - 2885154, 

94 = igA + 2AK-^gl. 

Eliminating gi and 54 yields the Chazy equation (|27p . From the expressions for 7' and 7" 
one can easily obtain two algebraic relations connecting gi and 54 to the 7. Hence one can 
obtain, by solving these algebraic equations, the general solution in this subcase. Since it 
is well known that the Halphen system is equivalent to the Chazy equation, one may also 
express the solution in terms of the uJi . One can also show, without first having to solve the 
equations, that the solution is semi-simple. 
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4.2.2. N = A ,a = ^ . In this case we take 

F = ^tlU + tlt2t3 - ^{t2t3)h{tA) + {4 9i{U) + t2ti 53 (t4) + 1% gi{U)] . 

The analysis of the WDVV equations gives the single system (there are other subcases which 
appear in the analysis, but these yield non-semi-simple solutions): 

(31) V = ^'-72i^5i' 

g'l = 279; -517' 

where X is a constant and 



93 = Kg^ , 54 



■ , 1 ■ 
9i - o5i7 



30 

On eliminating gi one obtains a third-order scalar equation 

... 3(7-277)2 .. 



7-7' 



-|- 677 — 67 . 



Here the independent variable has been rescaled, t = , and hence 7' = ^7 etc.. This falls 
within Bureau's class (see appendix) and its solutions are given in term of the Schwarzian 
triangle function y{t) = S[^, ^,0,t] . 



2 \ y y-lj 

ill f y ] 

2dt °^\yl(y-l)lj • 

Given this solution one may easily find the remaining functions: they all take the schematic 
form 



for various constants aj , 5j , Cj . 



4.3. N = 5 solutions. In this case there is again a 1-parameter family of Euler vector fields. 

In this case we have 3 pivot points which correspond to the monomials {(^2^4)^5 ^2(^3)^^47 (^s)^}- 
Again, without loss of generality we may assume decreasing degrees and hence o" < ^ . An- 
alyzing the Frobenius equation ()26p . we find the non-pivot terms summarized in table 2. 

Once our polynomial ansatz is made, WDVV reduces to systems of ordinary differential 
equations for the pivot and non-pivot functions. In all cases, except for a = 1/2, these 
equations typically fall into one of two types: 

1. second order Euler- type differential equations that may be integrated explicitly to 
yield a non-semi-simple solution to WDVV. 
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cr 


Non-pivot solutions {02,03,04): 


1 
2 


(4, 0,0), (3, 0,1), (1,0, 3), (0,0, 4) 
(3, 1,0), (2, 1,1), (1,1, 2), (0,1, 3) 
(2, 2,0), (0,2, 2) 
(1,3,0), (0,3,1); 


1 

3 


(3, 0,0), (1,0, 4), (0,0, 6) (0,2,3); 


1 

4 


(1,0,5) , (0,0,8) 
(2, 1,0), (1,1, 3), (0,1, 6) 
(0,2, 4), (0,3, 2); 


^ , (n > 5 , and odd) 


(l,0,n + l),(0,0,2n) (0,2,n); 


ra'("^2) 


(0,0,2n + l); 




(1, 0, 2n + 1) , (0, 0, 4n) , (1, 1, n + 1) , (0, 1, 3n) 
(0,2,2n),(0,3,n) 


lr'(^>4) 


(0,1, n) 


arbitrary 


none . 



Table 2. Non-pivot solutions of the Probenius equation ()26p for N = 5 . 



2. first order equations of the form 

(32) V - = 0, 

where 7 is a pivot term. 

In such cases the resulting Frobenius manifold is non-semi-simple and therefore a = 1/2 is 
the only value for which one may hope to find a semi-simple modular Frobenius manifold. 
Note that this is the most symmetric case - the prepotential contains 15 unknown functions 
(3 pivot and 12 non-pivot) and homogeneity means that the overdetermined system to be 
solved is far more complicated than in other examples. We note, though, two solutions which 
have appeared in the literature. 

Example 12. The following example (resettled slightly) was found in [lOj . 

F = ^tjh + h (t2U + \ti^ - \ (t2u + ^ti) 7(^5) 

+ ^ {4 - 2*2*1 - 2*2*1 + 3*1*4) hl{h) + ^ {4 - 2*4*1 - 2*2*4 + 3*i*l) /l2(*5) 
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where 

1 = ip -^hih2, 
h[ = 7/11 + hi , 
h'2 = 7/12 + hf . 
This system reduces to the Chazy equation for the function 7 . 

Example 13. The following example may be obtained by 'folding' the D^'^^ solution found by 
Satake [24J (see section\B^, reducing a six dimensional solution to a five dimensional solution. 

F = \tlh + \ti{tl + tl + t'^-^{tl + tl + tf,^^{t^) 

+ (4 + tt- 2tltl - 2tltl - 2tltl) (71 (ts) + ^(^2*3*1)52(^5) , 

where 

7' = ^7' - ^52' - 1285? , 

a'l = 191 - Y^gl + 1%? , 

92 = 192 - 325132 • 
This system too reduces to the Chazy equation for the function 7 . 

Note that the two systems in these examples are inequivalent - the corresponding quadratic 
algebras are non-isomorphic (see section [7] and jl7j ) . 

5. Modular almost-dual solutions 

In [8] Dubrovin introduced the notion of an almost-dual Frobenius manifold. Consider the 
vector field defined by the condition 

E-^oE = e. 

This is defined on M* = M\S , where S is the discriminant submanifold on which is 
undefined. With this field one may define a new 'dual' multiplication -k : TM* x TM* — )■ TM* 

by 

X i.Y = E-'^ o X oY , yX,Y£TM\ 

This new multiplication is clearly commutative and associative, with the Euler vector field 
being the unity field for the new multiplication. 

Furthermore, this new multiplication is compatible with the intersection form g on the 
Frobenius manifold, i.e. 

g{Xi<Y,Z)=g{X,Yi<Z), yX,Y,Z€TM*. 

Here g is defined by the equation 

g{X,Y) = rj{X oY,E-^), yX,YeTM* 

(and hence is well-defined on M* ). The intersection form has the important property that it 
is flat, and hence there exists a distinguished coordinate system {p} in which the components 
of the intersection form are constant. It turns out that there exists a dual prepotential F* 
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such that its third derivatives give the structure functions c*-^ for the dual multipHcation. 
More precisely [8]: 

Theorem 14. Given a Frobenius manifold M , there exists a function F* defined on M* such 
that: 

* - ( ^ ^ ^ 

dp'^dpidp^ 

Moreover, the pair {F*,g) satisfies the WDVV-equations in the flat coordinates {p} of the 
metric g . 

Thus starting with a Frobenius manifold one may construct an 'almost-dual' Frobenius 
manifold (for full details see [Bj) with prepotential F* : 

F ^ F* 



Combining this with the symmetry / leads to the following picture 




In |15) the induced involutive symmetry 

F* F* . 



was constructed. It turns out that the action I* is different for modular and non- modular 
Frobenius manifolds. 

Theorem 15. Let F define a Frobenius manifold and let F denote the induced manifold 
under the action of the symmetry I . Let F* and F* denote the corresponding almost dual 
structures. The I* , the induced symmetry act as: 

• Case /; d / 1 : 



9ab 

F* 



p_ 

il' 

El 



1,... ,N, 



where ti = \gahP°'P^ ■ 
Case II: d = 1 : 
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2 h ' 

^, i = 2,...,N-l. 



1 

9ab — 9ab ! 



Thus if one starts with a modular Frobenius manifold the almost-dual prepotential will 
also satisfy Proposition [6l 

In all known examples (essentially the Jacobi group orbit spaces for and Bn^2 for 

arbitrary N ) the (modular) almost dual prepotential takes the form 

(33) F* = \^u\ - ]^u{7., z) + ^ Kf ((a, z), r) 

where f{z,T) is constructed from the elliptic trilogarithm (see [27]) and has the expansion 
f{z,T) = -—^ I log z + log r/(r) 



{2m) [2 

+ (2^ § T^^^l^^^^"^'""'' ' 

where i?2n are Eisenstein series - and hence are modular forms not functions - and it follows 
from this that / has the required modularity property 

/ ( — , ) = —^f(z,T) ir— , (mod quadratic terms) . 

\T T J T'^ 4! 

The set of vectors 11 is constructed via a Landau-Ginzburg/Hurwitz space construction [18]. 
These satisfy the conditions 



^/iQ,(a,z)^ = 0, 
Y^ha{a,z)'^ = 3(z,z)^, 



aeiX 



and it then follows that the basic modularity property given in Proposition [U] is satisfied. 
Using these conditions the almost dual prepotentials may be written in a slightly modified 
version of Proposition [TUl 

1 1 2Tri 1 

F* = ^u\--uiz,z)-—E2iT){z,zf-^^Y,K{a,zflog{a.z) 



+ ^c„£'2„(t) <^ ^ha{a.z] 

n=2 loGil 



(2n+2) I 
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Note the single pivot-term, proportional to E2{t) . 

In |27) sufficient conditions were given on an arbitrary set of vectors il to ensure that a 
function of the form (|33p is a solution of the WDVV equations. However, determining further 
conditions which would ensure whether or not the solution is also the almost-dual potential 
of a modular Frobenius manifold (or even a polynomial modular Frobenius manifold) is a 
considerably harder problem. Another open problem is the construction of the almost-dual 
manifolds corresponding to the low-dimensional polynomial modular Frobenius manifolds 
constructed in the previous section. Such a construction should be deeply connected to the 
theory of Jacobi forms O [22] . 

6. Foldings and codimension 2-toroidal Lie algebras 

As noted in the introduction, a number of polynomial modular Frobenius manifolds have 
appeared in the literature before. In fact, some of the earliest solutions of the WDVV equa- 
tions fall into this class and were written down before the formal definition of a Frobenius 
manifold. The space of versal unfoldings of the simple elliptic singularities, denoted by dif- 
ferent authors as £"6,7,8 j -^eVs '-'^ -^8 j ^9 ) <^io > carries the natural structure of a Frobenius 
manifold. The prepotential and Euler vector field for the Eq singularity were derived by 
Verlinde and Warner [3Uj : 

F = ^tlt8 + h{t2t7 + hte+Uh)+t2t3Ufo{t8) + ^{tl+4+tl)fi{ts) 

+ (^2*3^6*7 + *2*4*5*7 + *3*4*5*6)/2(i8) + ^ (^2*5*6 + ^3*5*7 + i4*6i7)/3 (^s) 

+ {t2t3tl + t2Utl + t3Ut^j)U{ts) + ^(t^tf + tltj + tltDMts) 

+ ]:[t2t7i4 + 4) + ht6{4 + 4) + Ut5i4 + 4)]hits) 

D 

1, 2 2 2 xx/A'^/SS 22 2 2\ £ ( \ 

+ 2 1*2*5*6*7 + *3*5*6*7 + *3*5*6*7)/7(*8) + ^1*2^5^6 + *3*5*7 + *4*6*7)/8(*8) 

+ ^(*2*7 + *3*6 + *4*^)/9(*8) + ^(*i*i + *!*? + *^*?)/lo(*8) 

+ ^(*5*6*7 + *5*6*7 + *t*6*7)/ll(*8) + ^*i*6*?/l2(*8) + 7^(*5 + *6 + *7)/l3(*8) , 

and 

9? + 3 V 9^ + ' + ' 9Fj+3 V + ' 9^+* • 
The functions ji may all be expressed in terms of the Schwarzian triangle function i, i, t] . 
Note that /2 and are the pivot-terms. The prepotential for E-j was derived by [13]. The 
prepotential for E^, has not been written down - one would expect it to have over 200 terms. In 
principle this is just a computational exercise - one can find the multiplication in terms of the 
natural deformation variables in the Jacobi ring and then use the flat coordinates derived by 
Noumi and Yamada to derive the prepotential [16j. It is clear, though, that this prepotential 
will fall into the class of polynomial modular Frobenius manifolds considered here. It is 
interesting to note that the invariance of such solutions under an inversion symmetry was 
known already by Verlinde and Warner [30j , but the geometric origins of this symmetry were 
not fully understood. 
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The connection between the underlying singularity theory and the Frobenius manifold is 
subtle in these examples: one has a one parameter family of Frobenius manifolds associated 
with each elliptic singularity, and to construct a global theory one has to understand object 
such as the primitive form as a function of this parameter (the calculations in [30j and [16] 
were performed at one point in this family). For details of this global construction see the 
recent preprint of Milanov and Ruan [T^j. From the point of view of this paper it is hard see 
this one parameter family; the modular symmetry condition in Proposition [TD] may impose 
too strong conditions to see this family. One possibility might be to use the weaker version 
of modular symmetry used in Example HI 

Another polynomial modular Frobenius manifold was derived explicitly by Satake [24J and 
denoted D^^'^^ . Here the prepotential and Euler field are given by 

F = \tlh + \h{tl + tl + tl + tl) + ]^[tl+tl + tl + tfj^-i{t^) 

+ [4 + 4 + i\ + 4- 2(^2*1 + ilil + thl + t¥l + ^1*5 + ilil)] 9i{h) + (^2*3*4*5 ) 52(^6; 

id 1 / 3^ 4^ 5^ 

It was shown by Zuber [3T], in the case of polynomial Frobenius manifolds, that the non- 
simply-laced examples (namely Bn ,^4,^2, -^3,4 , hi^)) may be obtained by restricting the 
simply laced examples (namely An , Dn , -E'ej.s) to certain hyperplanes S = {i* = , i ^ /} for 
some subset / C {1 , . . . , A^} . In order to obtain a subalgebra one requires the condition 

k 



and 



^ = 0, e I ,yk <^ I . 



This process may be understood in terms of automorphisms of the corresponding Coxeter 
diagrams and their foldings (see, for example. Figure [1]). In such diagrams the arrow should 
be read as 'folds to' . 



Figure 1 . The folding of to F4 . 

Similarly the solutions found by Verlinde and Warner [30j and Satake [23] given above also 
admit such foldings, and the resulting solutions remain in the class of polynomial modular 
manifolds. 

The d'^'^^ submanifolds ar^: 

^2,a — {t2 = ts = ti} 
/ \ 
^ Si ^ {t4 = ^5} S3 ^ {t2 = H = t^= ts} 

\ ^ ^ 

5^2,fe — {^2 = ^3 5 ^4 = ^5} 



^Here the subscripts on E denote the codimension of the submanifold (i.e. the number of constraints on 
the original, ambient, manifold). The empty set denotes the absence of any constraints and so corresponds to 
the original manifold. As more foldings take place the number of constraints increase, and the dimension of 
the submanifold decreases. 



POLYNOMIAL MODULAR FROBENIUS MANIFOLDS 



19 



and the Eq ' submanifolds are: 

^ Si ^ {t2 = ts ,^6 = h} ^ {^2 = ^3 = ^4 ,t5 = ^6 = ^7} ■ 

The corresponding foldings are shown Figures [2] and [3l Note that the E^'^^ manifold does 
not admit a folding to a five dimensional manifold, in agreement with the results of section 
m These foldings may also be understood in terms of boundary unimodal singularities yai2j) 
and in terms of the theory of 2-toroidal Lie algebras \21\ |23] . More details on the meaning of 
these more complicated diagrams may be found in [2T]. However, there are subtleties in the 
associated invariant theory \26\ . 




Figure 3. The folding of . 

The folded prepotentials and the prepotentials found in section[5]do not, at first sight, agree: 
the underlying triangle functions and the forms of the functions that appear are different. 
However, the two different triangle functions are related by certain nonlinear transformations 
which have their origin in the nonlinear Goursat identities for hypergeometric functions 
For example, the relation 

leads to the relation x = between the Schwarzian functions x = ,0,0, (— 4) et] and 

X = , ,0 ,t] . 



7. Modular dynamical systems 

The form, and modular invariance, of the equations derived in section [l] suggests that 
one should develop a theory of modular dynamical systems. Suppose that one has a set of 
quasi-modular functions (we assume no analytic properties such as being holomorphic in the 
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upper-half plane) 7(t) and ^^(t) , a G W for some indexing set W such thalS: 
7 ^ — ^ = T^7(r) + ar , a / constant , 

Such a Qn is said to have weight n . Define next a Rankin-type derivative 

^ dgn 

^{9n) = 9n , 

ctr a 

^ 7 = 7 • 

ar a 

It is easy to check that -D(7) has weight 4 and D[gn) has weight n + 2 and that D{gngm) = 

gnDgm + gmDgn ■ 

Definition 16. Let W he some (finite) indexing set. A modular dynamical system takes the 
form 

D{l) = q{g), 

D{ga) = Pais), a£W. 

where the polynomials q and are constrained so that the resulting system is invariant under 
the transformation r ^ — ^ . 

The interplay between the weights and the polynomials places very strong restriction on 
the form of these equations. Note that if the polynomials depended on 7 as well, then the 
modularity property would be lost. Thus 7 only enters through its appearance in the Rankin 
derivatives. 

Example 17. Consider the 3-component system with W = {n , m} . Clearly if m > 5 and 
n > 5 then = and one obtains a degenerate system (we say an n-component system 
is degenerate if the differential equation for 7 is of order less than n ) and a similar analysis 
shows that the systems {4, m > 7},{3, m > 6},{2,m > 5},{l,m > 5} are also degenerate. 
Thus there can only be a finite number of non- degenerate systems. Certain other cases are 
degenerate, and the final list of non- degenerate cases is: 

W = {6, 4} , {4, 3} , {4, 2} , {4, 1} , {3, 1} , {2, 2} , {2, 1} , {1, 1} . 

Consider now the system W = {6,4}. Then the modular dynamical system must take the 
form: 

D{-f) = ag4 , 
D{gA) = bge , 
D{g6) = cgl 

for constants a ,b ,c. A famous example in this class is the Ramanujan differential equations 
for the Eisenstein series {E2 ,£'4 ,-^6} • Eliminating g^ and ge gives the third order system 

D^j + k{Djf = 



■^The value of the non-zero constant a is not fixed here - this is to facihtate comparison with the work of 
various sets of authors who use different values in their work. 



POLYNOMIAL MODULAR FROBENIUS MANIFOLDS 21 

where k = — ^ (if a = the system is degenerate). This system again falls into the class of 
generalized Chazy equations considered in the appendix and its solutions can be expressed in 
terms of Schwarzian triangle functions. 

The differential equations obtained in section [5] may be written as modular dynamical 
systems: 

= 4 , fj = ^ Case I: The system (|28p may be written as: 

D{h) = -4g\ D{g)=h, D{j) = g\ 
Here the weights of g{= gi) and h are 2 and 4 respectively. 

= 4 , 0" = ^ Case H: The system ()30p may be written as: 

D{g) = {2AK)h^ , D{h) = {2AK-^)g^ , D{^) = -288gh . 
Here the weights of g (= gi) and h (= (74) are both 2. 

A'^ = 4 , fj = I The system ()3ip may be written as: 

D{g3) = ^g!, D{gi) = g3, D{j) = -gf . 

Here the weights of gi and 53 are 1 and 3 respectively. 

There is considerable scope for reducing such modular dynamical systems to canonical form 
by redefining, if the weights allow, the functions 7 and gn • 

Example 18. Consider the system with W = {1 , 3} . The most general system is 

Dgi = agf+bgs, 
(34) 1)7 = cgf + dgigs, 

Dgs = egf + fgfg3 

where a ,b , . . . , f are arbitrary constants. Within this system one may make the following 
transformations 

93 ^ 93 + cigf , 
7 7 + C29I 

and with a careful analysis of the various (ignoring degenerate) cases one obtains the following 
canonical forms: 

6 = 0: Dgi = 0, 

D-f = gigs , 

Dgs = agl + /Sgjgs , 

6 / : Dgi = gs, 

D7 = -9f, 

Dgs = agl + /Sgfgs , 

where a and /3 are constants. 

On eliminating gi and g^ one obtains a third order modular ordinary differential equation 
for 7 from whose solution one can obtain gi and 53 by quadrature. The third order system 
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will not, in general, fall into Bureau's class, and the requirement that it is in this class will 
place restrictions on the parameters. 
For the W = {1 , 3} system above: 

• If b = then it falls into Bureau 's class if and only if a = , 13 ^ ora^0,/3 = 
(if a = l3 = the system is degenerate). The corresponding systems are, respectively 
{D-f){D^-f) = {D'^-ff and D^-y = . 

• If b ^ then it falls into Bureau 's class if and only if f3 = and the corresponding 
system is 

{D^){D-'^) = ^{D'jf-Aa{D^f. 

The solution is given in terms of the Schwarzian triangle function S'[^, ^,p,t] where 
= I — . Note that if one requires the global condition that the triangle functions 
form a arithmetic triangle group then this forces p = and hence a = | which is 
precisely the system (f3T]l . 

An alternative method to solve such modular dynamical system is to use a nonlinear change 
of dependent variable. Again, for the {1,3} system one may define: 

91 

With this the system ()34p becomes 

X = X^ + cY^ + dYZ , 

Y = 2XY + 2aY^ + 2bYZ , 

Z = 2XZ + eY"^ + (/ - a)YZ - bZ"^ . 

This falls into the class of (three-component) quadratic dynamical systems considered by 
Ohyama [17| : 

In fact, all the above examples (except possibly the {1, l}-system) may be written, via ap- 
propriate nonlinear changes of variable, in such a quadratic form. The method of solution 
depends on the commutative algebra defined by 

Xi o Xj — ^ ^ ^ij^k • 
k 

In particular the modular invariance of the dynamical systems means that the corresponding 
algebra carries a unit. Once the algebra is identified (a complete classification of three dimen- 
sional algebras with unity was carried out in [T7]) the solution is given in terms of solutions 
to second order Fuchsian ordinary differential equations. In all examples appearing in this 
paper this turns out to be the hyper geometric equation. 

We note also that in all existing examples the underlying modular dynamical system is of 
rank three. As dimension increases one might expect to find higher order dynamical systems 
underlying the WDVV equations, but in all the examples so far found one just obtains third- 
order systems. 
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8. Comments 

As mentioned in the introduction, some of the earhest examples of solutions of the WDVV 
equation [3^ [T^ were actually examples of polynomial modular Frobenius manifolds. Re- 
cently these have attracted renewed attention due to their appearance in orbifold quantum 
cohomology, i.e. the Gromov-Witten theory of P^^ (the complex projective space with 
ri-orbifold points). In particular one has an isomorphism (mirror symmetry) between the 
orbifold picture and the singularity theory picture |T21[2S], e.g. 



The genus-1 contribution to the free energy may also be calculated for these examples [28j . 
In fact, the modularity property inherent in the genus-0 free energy also manifests itself at 
genus-1 [22] • Recently Dingdian and Zhang [B] have conjectured how higher- genus correction 
should transform under inversion symmetry. Ultimately one would wish to understand how 
the full dispersive hierarchy transforms under this modular symmetry: currently only the 
genus-0 (or dispersionless) integrable systems transformation properties are known [6l I15j. 
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Appendix 

In this appendix we summarize the results of Bureau who constructed the general 
solution to the third order differential equation 

(35) u=a- ^ + I3uu + '-^u'^ + 5{u - u^Y 

u — 

in terms of Schwarzian triangle functions. This will be referred to as Bureau's class. Here the 
constants a , /3 , 7 , (5 are constrained by the relations 

16a -7 = 18, 2/3 + 7 = 6. 

The solution is given in terms of a function y , namely 

^ = M ? - Q{y)y 
2 [y 



where 



y y - i 

and y satisfies the Schwarzian differential equation 

y 3[j/]^_ l[ l-n^ 1-m^ -m^ -n^) y 2 

^ ^ y 2\yj - 2I y2 +(i_y2) }y 
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The relationship between a , /3 , 7 , (5 and m ,n ,p are given by the requirement that Lj = , i = 
1,2,3 where 

Li = (1 - 2n) [1 - 3n - a(l - 2n)] , 

L2 = (1 - 2m) [1 - 3m - a(l - 2m)] , 

L3 = (1 - 2m)(2 - 3n) + (1 - 2n)(2 - 3m) - 2a(l - 2m)(l - 2n) 

+ ^(7 + J - 6) [(1 - m - n)2 - p2] . 

As is well known, the solutions of the Schwarzian equation (j36p are given in terms of 
the Schwarzian triangle function y = S[m,n,p,t] which may be constructed as the inverse 
function of the ratio of two solutions to the hypergeometric equation with parameters a ,b ,c 
(i.e. one solution is 2Fi[a,b, c, z]) where 

= (1 - c)2 , = {a + b- c)^ , / = (a - b)^ . 

Thus starting with a linear equation - the hypergeometric equation - one may, by taking ratios 
and inverse functions, construct solutions to Bureau's equation (|35p . In terms of the Rankin 
derivatives ()35p may be written as 

{Du){D\) = a{D^uf + {6 + 16a - 24){Duf 

with /3 = 12 — 8a , 7 = 16a — 18 . Note that with these formulae: 



(1 — m — n)^ — p 



y{y - 1) 

In fact. Bureau's method is considerably simplified on the use of the Rankin derivative. It 
is easy to show that 

where hn+i{y) = ^^^^ + (n+l)g(y)/i„(y) . Thus Bureau's equation ([55]) becomes the algebraic 
equation hi{y)h2,{y) = a/i2(y)^ + fihi{y)^ which should hold for all values of y . Equating 
coefficients gives the three equations Lj = . There is considerable scope for generalizing 
this method by replacing the hypergeometric equation with more general Fuchsian or even 
non-Fuchsian systems: again one just obtains algebraic conditions. 

When a = one obtains Chazy's equation class XII with solutions given in terms of the 

function S[\,\,p,t\ where = ^ (5^4) ™^ Chazy equation itself occurs when, in 

addition, 5 = ^ with solution given in terms of the function |, 0, t] . 

A systematic procedure to solve such modular differential equations has been derived in 
[5] . It is interesting to note, though, that all the examples of modular Frobenius manifolds 
constructed here come from solutions of Bureau's equation and not more general equations 
(or, more notably, higher-order modular equations). 
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